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Abstract. We investigate the effect of varying boundary conditions on 
the renormalization group flow in a recently developed noncommutative 
geometry model of particle physics and cosmology. We first show that 
there is a sensitive dependence on the initial conditions at unification, 
so that, varying a parameter even slightly can be shown to have drastic 
effects on the running of the model parameters. We compare the run- 
ning in the case of the default and the maximal mixing conditions at 
unification. We then exhibit explicitly a particular choice of initial con- 
ditions at the unification scale, in the form of modified maximal mixing 
conditions , which have the property that they satisfy all the geometric 
constraints imposed by the noncommutative geometry of the model at 
unification, and at the same time, after running them down to lower 
energies with the renormalization group flow, they still agree in order of 
magnitude with the predictions at the electroweak scale. 
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1. Introduction 




The recent work [26] developed cosmological models of the very early 

universe based on the particle physics model of |13j derived from noncom- 
mutative geometry, via the formalism of spectral triple and the spectral 
action functional. Other results on cosmological aspects of noncommutative 
geometry models of particle physics include P]. P], [?7|- PHj- PI- 

In the particle physics model of |1] , [T^] , jlSj , the Lagrangian is obtained by 
computing the asymptotic expansion at high energy of the spectral action 
functional [TT] on a noncommutative space which is the product of an ordi- 
nary (commutative) spacetime manifold and extra dimensions given in the 
form of a noncommutative space which is metrically zero-dimensional, but 
K-theoretically six dimensional. The choice of the noncommutative space 
determines the particle physics content of the model and the gauge sym- 
metries. The masses and mixing angles arise geometrically as coordinates 
on the moduli space of Dirac operators of the spectral triple describing the 
extra dimensions. In the case of the model developed in [13], the particle 
physics content is the same as in the i/MSM, namely, in addition to the par- 
ticles of the Minimal Standard Model, one has right handed neutrinos with 
Majorana mass terms. However, the model is significantly different from 
z^MSM when it comes to the properties of the action functional. In fact, as 
proved in [13] . the asymptotic expansion of the spectral action contains the 
full Standard Model Lagrangian, with the additional Majorana terms for 
the right handed neutrinos. One has unification of the coupling constants of 
the three forces, hence the model has a preferred energy scale at unification. 
The asymptotic expansion of the spectral action also contains gravitational 
terms, which are the most interesting part from the point of view of applica- 
tions to cosmological models. These terms contain an Einstein-Hilbert term, 
a cosmological term, a conformal gravity term, a non-dynamical topological 
term, and a conformal coupling of the Higgs field to gravity. 



In the approach to cosmological models developed in [26], one uses the 

fact that, at the unification scale, in the terms one obtains in the asymptotic 
expansion of the spectral action, the usual gravitational and cosmological 
constants are replaced by effective constants, whose expression at unification 
depends upon the Yukawa parameters of the particle physics content of the 
model. In [26j . this is used to derive an early universe model in which one 
allows these effective gravitational and cosmological constant determined 
by the boundary conditions given by the asymptotic expansion of the spec- 
tral action, to run with the RGE fiow of the associated particle physics 
model, according to the running of the Yukawa parameters and Majorana 
mass terms. This allows for a much more serious variability of the effective 
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gravitational and cosmological constant in between the unification and the 
electroweak epochs of the very early universe (and in particular during the 
inflationary epoch) than is usually considered in other gravity models. In 
, this type of running leads to several consequences on early universe cos- 
mology, from mechanisms for inflation to effects on the gravitational waves 
and the evaporation law for primordial black holes. 

For the purpose of the present paper, the specific issues of the running of 
the gravitational parameters and of the resulting interpretations within the 
model are not directly relevant, since the results we give here are specifically 
about the running with the RGE flow of those expressions of Yukawa pa- 
rameters and Majorana masses, which enter the value at unification of the 
asymptotic expansion of the spectral action. 

The analysis performed in [26] depends on the choice of initial boundary 
conditions at unification for the renormalization group flow. The results of 
[26j are obtained using the default boundary conditions of [Ij. However, as 
we show in the present paper, one obtains signiflcantly different behaviors 
of the coefficients of the asymptotic expansion of the spectral action by 
changing boundary conditions. This implies that there will be the possibility 
of drawing interesting exclusion curves in the space of all possible boundary 
conditions, on the basis of comparing the model with cosmological data, 
for example through the predictions for the tensor-to-scalar ratio and the 
spectral index derived in [26]. 

For the purpose of the present paper, we first show how one obtains 
significantly different curves for the running of the parameters in the action 
functional with different choices of the boundary conditions. This shows, as 
one would have expected, a sensitive dependence on the initial conditions at 
unification, which means that a fine-tuning problem arises within the model, 
in the choice of the data at unification. 

The main result of the paper is then to exhibit a specific choice of bound- 
ary conditions, which we denote modified maximal mixing conditions, which 
differ from the default one of p[], and which have the desired properties. 
Namely, we show that all the geometric constraints on the data at unifica- 
tion derived in \13} are satisfied by our choice of boundary conditions. We 
also show that, when running the RGE fiow with those boundary condi- 
tions, one obtains values in the low energy limit that are still compatible in 
order of magnitude with the physical predictions and observed values at low 
energy, as in the case of the default conditions of [1] . 

An important aspect of these models is understanding how much nonper- 
turbative effects in the spectral action may affect the low energy behavior of 
the model, since that is the main obstacle to extending to the more recent 
universe the cosmological models of |26j . Our estimates of the low energy 
behavior when matching geometric boundary conditions at unification may 
also provide some indirect evidence for the magnitude of such effects. 

Recent results of [27J show that, at least in the case of sufficiently symmet- 
ric geometries, the spectral action can be fully computed non-perturbatively. 
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using the technique of [12], and the non-perturbative effects are hmited to 
the shape of the inflation potential. 

2. The spectral action and the renormalization group flow 

In noncommutative geometry one models the analog of a Riemannian 
manifold through the notion of a spectral triple, consisting of data {A, %, V) 
of an involutive algebra, a Hilbert space representation, and a Dirac opera- 
tor, which has the compatibility condition of having bounded commutators 
with elements of the algebra. Additional structure, in the form of grading 7 
and real involution J with compatibility conditions with the data {A, Ti, T>) 
are also introduced. In the particle physics context, 7 corresponds to the 
two chiralities of fermions and J to the involution that exchanges particles 
and antiparticles. See [13] for a more detailed account of the underlying 
mathematical structure, which we do not recall here. The action functional 
considered in noncommutative geometry models for particle physics is based 
on the spectral action jllj for the Dirac operator of a spectral triple, with 
additional fermionic terms. In the model of [13] this takes the form 



(2.1) TT{f{DA/^)) + \{JlDAi). 

Here Da = D + A + e' J AJ^"^ is the Dirac operator with inner fluctuations 
given by the gauge potentials of the form A = = Ylk'^k[D,bk], for ele- 
ments ak,bk G A. The e' is just a function of n mod 8 that gives -1 for n 
congruent to 1 mod 4 and 1 for all other values of n. The fermionic term 
{J(,,DaO should be seen as a pairing of classical fields G Ti'^ = {S, £ 
T-L\jS, = C}) viewed as Grassman variables. For the purpose of cosmological 
applications, the most important part of this action functional is the one 
that comes from the asymptote expansion at high energy A of the spectral 
action Tr(/(L'^/A)), since this contains the gravitational terms and their 
coupling to matter. 

2.1. The asymptotic form of the spectral action. The asymptotic ex- 
pansion of the spectral action is obtained in the form (see jllj . |13| ) 

(2.2) Tr(/(I)/A))~ Yl /fcA'/l^r' + /(0)Ci^(0)+o(l), 

A;eDimSp+ 

where fk = f{v)v^~^dv are the momenta of the function / and the 
noncommutative integration is defined in terms of residues of zeta functions 

(2.3) Ca,D{s) = Ma\D\~n- 



The sum in (2.2) is over points in the dimension spectrum of the spectral 
triple, which is a refined notion of dimension for noncommutative spaces, 
consisting of the set of poles of the zeta functions ( |2.3| ). More explicitly, 
as proved in [13], when applied to a noncommutative space of the form 
X X F, with X an ordinary 4-dimensional (Euclidean) spacetime and F 
the noncommutative space whose algebra of coordinates is C ® H© M3(C), 
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with H the algebra of quaternions, the expansion (2.2) of Tr(/(Z)yi/A)) gives 
terms of the form 

+ aoj C^.p^C^^^'^^d^x + ro j R*R*^d^x 

(2.4) + ^j\DH\^^d''x-filj\H\''^d''x 

- (oj R\H\^^d^x + Xo I m^^gd^x 

+ \ j (G'pu G''"' + F'^-" + B,, ^ d^x. 



The coefficients of these terms are functions 

96/2 A2 - foe 



(2.5) 
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These depend upon the three parameters fo, /2, fi, where fo = /(O) and 
for /c > 

fk = / /(^^)i^'-'d^, 
Jo 

where fo depends upon the common value of the coupling constants at uni- 
fication energy and /2 and /4 are free parameters of the model. The expres- 



sions (2.5) also depend upon the energy scale A and the running of these 



parameters is the main topic of our present investigation. In addition to the 



explicit dependence on A of the coefficients (2.5) there is also an additional 



and very interesting dependence on A through the coefficients a, b, c, d and 
e. These are functions of the Yukawa parameters and Majorana masses of 
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the particle physics content of the model, in the form 



(2.6) 



I) 



a 



b 



e 



c 



Tr((yjy,)2 + {Yhef + 3(yjy.)2 + siY^Y^f) 

Tr(MMt) 

Tr((MMt)2) 

Tr(MMtyJyj,). 



2.2. Renormalization group flow. The particle physics models based on 
the spectral action functional of noncommutative geometry as in |13] 
are (at present) entirely a classical theory. In particular, this means that 
whenever physical predictions are derived in these models using renormal- 
ization group techniques to lower the energy scale from unification, where 
the model naturally lives, to ordinary energies, one uses beta functions and 
renormalization group equations that are imported from the ordinary QFT 
of the specific particle physics Lagrangian that is obtained from the as- 
ymptotic expansion of the spectral action. This is a delicate issue, since in 
fact the asymptotic expansion includes both matter and gravitational terms. 
The beta functions and RGE fiow adopted here (as in |26j ) is the one for 
the extension of the Minimal Standard Model that includes right handed 
neutrinos with Majorana mass terms, while the gravitational effects are not 
taken into account in the form of RGE. This is an approximation, since the 
non-minimal coupling of the Higgs to gravity in the model means that one no 
longer has a clear separation between the particle and gravitational sectors. 
Consequences of modified RGE flows coming from non-minimal couplings 
to the Higgs can be found for instance in [7|, [8], and in [32], while effects 
from gravity terms are considered in jl8j. For the Minimal Standard Model, 
there is an extensive literature on the form of the beta functions and the 
RGE flow, see for instance [23] and references therein. In the case of the 
noncommutative geometry model of particle physics of [H], which did not 
yet include right handed neutrinos and Majorana mass terms, predictions 
of the Higgs mass were obtained based on using the RGE of the Minimal 
Standard Model. 

The RGE analysis of the model of [T3] considered in [26], which we also 
work with in this paper, differs from the usual RGE analysis of the Standard 
Model in the following ways: 

• Instead of the RGE of the Minimal Standard Model, one considers 
the equations for the extension with right handed neutrinos and 
Majorana masses, as in [T]. As in [T] these are treated by considering 
different effective field theories in between the different see-saw scales 



(see also [2], [3]). 
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• We vary the initial conditions at unification, by imposing the geo- 
metric constraints derived in [13] and at the same time requiring that 
the low energy values remain close to the expected physical values. 

The specific information on the NCG model of [13] enters here in two ways: 
first in selecting the appropriate matter content of the model (the presence 
of the extra right handed neutrinos with Majorana mass terms in addition 
to the usual Standard Model), hence the use of the RGE flow of [T], and 
also in the geometric constraints imposed on the boundary conditions at 
unification. 

We use, as in [2^] the renormalization group equations for the Standard 
Model with right handed neutrinos and Majorana mass terms of [IJ. The nu- 
merical results described here are obtained with a Mathematica code based 
on the REAP program of [1] adapted to our model by the first author. 

We recall here that the RGE for this particle physics model is given (at 
one loop) by the beta functions [Ij 

IGvr^ Pg^ = bigf with (&5C/(3)> &st/(2), ^{/(i)) = 15^ 

16^2 fs^^ = y,(^y]y, - ^yjy, + a-^gl- ^gj - Sgl) 
levr^ = y,(^yjy. - ^y,ty, + a-^g!- ^^i) 
levr^ = y,(^y,ty, - ^yjy. + a - ^5? - lal) 

167r2 /3m = y^yjM + M(Y^Y^f 
IGvr^ = - 3A(3gi + ^gf) + 3gl + ^{^gf + glf + 4Ao - 8b. 

Notice that we use here the normalization of the coupling constants used in 
[1], which is different from the one of |13j . 

In particular, as in [1], we solve numerically these equations using different 
effective field theories in the intervals of energies between the three see-saw 
scales, with matching boundary conditions. Namely, starting from assigned 
boundary conditions at unification, one runs the RGE flow down until the 
first see-saw scale (the top eigenvalue of the Majorana mass matrix M. Then 
one integrates out the higher modes by introducing a first effective theory 

(3) 

with Yukawa parameters y) obtained by removing the last row of Yy in the 
basis where M is diagonal and with Majorana mass matrix M(3) obtained 
by removing tha last row and column. One then restarts the RGE flow for 
these new variables with matching boundary conditions at the top see-saw 
scale, until the second see-saw scale, and so on. One has in this way effective 
field theories (yj''^ M(^)), /c = 3, 2, 1. 

We study the effect on this RGE flow of changing boundary conditions 
at unification scale, and we then derive consequences for the running of the 



8 



DANIEL KOLODRUBETZ AND MATILDE MARCOLLI 



coefficients o, b, c, d, e of (2.6). In the next section we show, as could have 
been expected, that the running is highly sensitive to the choice of the intial 
conditions at unification. This shows that there is an important fine-tuning 
issue in the model related to the assigned values at unification. We then 
present in the following section a specific choice of boundary conditions that 
meets all the geometric constraints on the model and that produces realistic 
values at low energies. 



2.3. A remark on gravitational and Yukaw^a parameters in the 
NCG models. This subsection is not directly relevant to the main result 
of the paper, which is simply a statement about the running of the parame- 
ters 0, b, c, D, e of (2.6), subject to different choices of boundary conditions 
at unification, with particular attention to those dictated by the geometric 
constraints imposed by the model of [13] at unification. However, we include 
it here to discuss briefly and compare different existing points of view on the 
role of the parameters (2.6) in the coefficients (2.5) of the spectral action 
expansion. 

In the NCG model of [13], the relation (2.5) between the coefficients of the 
asymptotic expansion of the spectral action and the Yukawa coupling and 
Majorana mass terms of the particle physics sector holds only at uniflcation 
energy. In particular, the dependence of the effective gravitational and ef- 
fective cosmological constants upon the parameters a, b, c, t), e of (2.6) only 
sets the boundary conditions at unification. In [13] (see also the exposition 
in Chapter 1 of [l6]), consequently, the running of the gravitational terms 
of the model is deduced from the usual approach as in |18j, see also [19J, 
by which one obtains only a very moderate (or essential lack of) running of 
the gravitational parameters. The running of the particle physics sector is 
then ruled, in the NCG models, only by the RGE flow of the matter La- 
grangian, neglecting gravitational effects (with the caveat mentioned above 
on the non-minimal coupling with the Higgs). 

However, there are cosmological models that include the possibility of a 
much more drastic variability of the gravitational parameters in the very 
early universe, including in particular the inflationary epoch. Scenarios 
with variable gravitational constant had been considered early on in Jordan- 
Brans-Dicke gravity, where the variability happens through the non-minimal 
coupling of gravity to a scalar field, and more recently within other modified 
gravity models, and in terms of RGE running [21], as well as in the context 
of primordial black holes with gravitational memory (see for instance [5J, 
or the recent [lOj and references therein). Similarly, a variable cosmological 
constant plays a role in various models (see for example [6], [20], [22], [33j). 

In [2n], therefore, a different viewpoint on the effective gravitational and 
cosmological constant in the asymptotic expansion of the spectral action 
in the NCG models is proposed, and a possible early universe model is 
investigated, which only covers the epochs in between the unification and 
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the electroweak eras, a period which is expected to include the inflation- 
ary epoch. It is shown that, if one considers an effective action where the 
gravitational and cosmological constant are allowed to run according to the 



RGE flow of the coefficients (2.6) through the expressions (2.5) and with 
the assigned boundary conditions at unification, then one recovers many 
of the scenarios predicted by other models with variable gravitational and 
cosmological constant, as [10], [1^, [221, and several different mechanism for 
inflation, with predictions about parameters such as the spectral index and 
tensor-to-scalar ratio. 

Other recent cosmological applications of [T^, such as those in |9j, [27] . 
[28], [29], [30], [31], follow the more conventional point of view on the as- 
ymptotic expansion of the spectral action and the form of the coefficients 



(2.5). These different viewpoints do not directly affect in any way the re- 
sults of the present paper, and we only mention them here for the reader's 
information. 



3. Effects of changing boundary conditions 

The REAP program from [Ij allows the user to adjust the boundary 
conditions. These changes are generally made at Aunif, taken here to be 
2 X 10^^ GeV. As we understand that only fine tuned initial conditions for 
the universe allowed its current form, we expect the boundary conditions at 
unification energy to drastically effect the development of our model param- 
eters. We show here, as an example, the different running of the coefficients 



a, b, c, d, e of (2.6) for the default boundary conditions and for the maximal 
mixing case. We also show explicitly the changing behavior of the running 
of one of these coefficients when one of the parameters varies at unification, 
in order to illustrate the significant dependence on the initial conditions. 



3.1. The default boundary conditions. The boundary conditions at uni- 
fication used in [26] are the default boundary conditions of These have 
the following values. 





YuiKnif) = I 0.00156368 

0.482902 
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For Yd{Aunif) = {Vij) they have 

yii = 0.0000482105 - 3.382 x lO-^^i 

yi2 = 0.000104035 + 2.55017 x 10"^^ 

yi3 = 0.0000556766 + 6.72508 x lO^^i 

y2i = 0.000104035 - 2.55017 x lO'^i 

= 0.000509279 + 3.38205 x IQ'^^i 

= 0.00066992 - 4.91159 x lO-^i 

= 0.000048644 - 5.87562 x lO^^i 

^32 = 0.000585302 + 4.29122 x 10"^^ 

yss = 0.0159991 - 4.21364 x lO-^^i 




Figure 1. The running coefficients with default boundary conditions near 

top see-saw scale 
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3.2. Maximal mixing example. To look at the maximal mixing case, 
we simply change Ynu at unification energy. With maximal mixing, our 
parameters will take these values. 



C = exp{2-Ki/3) 




UpMNsi-^unif) 



From the available estimates of the neutrino masses, we get the diagonal 
mass matrix 

12.2 X 10-^ 



3(ti) 






170 X IQ-^ 

15.5 X 10-3 



Finally, 



Yu = Ul,j^,j^g6(^^i)UpMNS 



Using this form for 1^ and the default boundary conditions on all the other 
parameters, we can look at the running coefficients. From the figures below, 
we see that there are vast differences in the development of the parameters 
with this boundary condition change. 




Figure 2. The running coefficients with maximal mixing boundary 
conditions near top see-saw scale 
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3.3. Running coefficients with changing boundary conditions. It is 

possible to get even more interesting behavior by using less standard bound- 
ary conditions. By changing just one parameter we can examine how it 
affects the flow of our running parameters. A specific example is the 1^ 
matrix. Using our standard boundary conditions, this matrix is diagonal at 
unification energy. We can adjust each of these elements on the diagonal, 
which correspond to our neutrino masses, to affect our flow. Using anima- 
tion functions in Mathematica, it is possible to get a nearly continuous idea 
of how the flow of our parameters develops with our boundary conditions. 
The figures below illustrate such a development discretely. 




Figure 2. Coefficient c at the upper see-saw scale with the first term of 
as 0.5, 1.0, and 1.5 respectively 

In these diagrams, we notice the transition changing as the upper neutrino 
mass varies. The sharp transition at the upper see-saw scale comes from the 
program integrating out the heavy neutrino at this scale. The second plot 
shows the behavior we expect from the standard conditions. In the first 
plot we can see the upper and middle transitions are much closer together 
than in our second plot. The final plot shows the transition at a much higher 
energy, corresponding to the higher neutrino mass. From these and other 
such plots, we learn how the running develops independently by changing 
different parameters. Of course, changing multiple parameters complicates 
this development and is dealt with in more detail when matching specific 
boundary conditions. 



4. Geometric constraints at unification 

There are some constraints on the boundary conditions at unification that 
are imposed by the underlying geometry of the model. These are derived in 
|13| . see also the discussion in §1 of |16] . We recall them here. Not all of 
these constraints are satisfied by the default boundary conditions of pQ, so 
a first improvement on the model of [26] is to identify choices of boundary 
conditions that satisfy these constraints, and then, among them, eliminate 
those that produce non-physical predictions. 



BOUNDARY CONDITIONS OF RGB FLOW IN NCG COSMOLOGY 



13 



We show here how to obtain a choice of boundary conditions that satisfy 
all the constraints by modifying the maximal mixing conditions. 

4.1. Constraint on A. A first constraint imposed by the geometry is on 
the value of the Higgs self-coupling A at unification. This satisfies 

Looking at our maximal mixing boundary conditions we can calculate 
that X{Aunif) = 2.989. By setting it to this value at unification energy in 
our flow we can ensure that this requirement is met. 

4.2. The a parameter and the Higgs vacuum. The model of |13j also 
relates the parameter a to the Higgs vacuum through the relation 



(4.2) 

TT g 

where g is the common value of the coupling constants at unification and 
M]y is the W^-boson mass. As M\y is directly proportional to ^/a, this 
condition is a statement of the equality of /o and the coupling constants at 
unification energy. 

4.3. Constraint on c. The see-saw mechanism is implemented in [T3] geo- 
metrically, through the fact that the restriction of the Dirac operator D{Y) 
to the subspace oiT-Lp spanned by vr, vl, Ur, Ul is of the form 








Ml 


Ml 
























Mr 








Ml 




v 








M^ 





) 



(4.3) 



where M^, is the neutrino mass matrix, see Lemma 1.225 of p!6]. This 
imposes a constraint at unification on the coefficient c, of the form 

(4.4) < C(A„„.,) < ""J^. 

Jo JO 

By setting our Majorana mass matrix to 10 times its default value, the 
inequality can be matched. In this particular case, the /2 that is used is 
in the range given in [26j. /o is calculated from the coupling constants at 
unification energy. 

4.4. The mass relation at unification. Another prediction which is spe- 
cific to the model of |I3j is a quadratic relation between the masses at 
unification scale, of the form 

(4.5) Yl + + + 3m^)|A=A„,„, = 8M^|a=a_,, 

generations 

where nii,, rrie, m„, and are the masses of the leptons and quarks, that 
is, the eigenvectors of the matrices S^i, 6^i, 5^^ and 6^^, respectively, and 
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Mw is the W-boson mass. We use the fact that Mw is given as a function 
of the model parameters by 

(4.6) y^ = Mw. 



2V2 



So, our equation becomes 



(4.7) ("^' + "^e + + 3m^)|A=A_^ = a|A=A. 

generations 

In our maximal mixing boundary conditions, we get 

(4.8) Yl + "^e + ^rnl + 3mi)|A=A„„,, = 0.6698 = a|A=A„„,,. 

generations 

This value of a, when converted to conventional units, gives a value of AIw 
of 72 GeV. The expected value on Mw is around 80 GeV so these boundary 
conditions are believable. 

4.5. Modified maximal mixing. Thus, the conclusion of this analysis is 
that we obtain a choice of boundary conditions that satisfies all the geomet- 
ric constraints of the geometric model at unification by using our maximal 
mixing boundary conditions as described in the previous section, but with a 
modified Majorana mass matrix and Higgs parameter, as explained here. We 
refer to the resulting boundary conditions as the modified maximal mixing 
conditions. 

We then need to check that, when we run the RGE flow with these initial 
conditions at unification, we obtain values at low energies that are compati- 
ble, within order of magnitude, with the required physical values. We discus 
this in the next section. 



5. Low ENERGY PHYSICAL CONSTRAINTS 

At the electroweak scale, physical data impose other boundary conditions 
on some of the Yukawa matrices. Finding the unification scale conditions 
that can also match these lower energy requirements is crucial to the theory. 
. We look at the conditions that are expected from physical data and com- 
pare to the results from the running of the model parameters. We show that 
our modified maximal mixing boundary conditions also satisfy the required 
constraints at low energy. 

5.1. Boundary conditions at the electroweak scale. Current predic- 
tions at the electroweak scale tell us that the CKM matrix at Aq can be 
taken to be of the form 

/ 0.97419 0.2257 0.00359 \ 
UcKMiAo) = 0.2256 0.97334 0.0415 

\ 0.00874 0.0407 0.999133 / 
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Combined with 

■ 0.00475 

246 



%3)(Ao) = ;^l 1 



4.25 

we get that the Yukawa parameters for the quarks are given by 
Yd = Uckm{^o)S(13){Ao)Uckm{MY 

and 

Similarly, for the matrix of charged leptons, the known values and low 
energy are 

' 0.000511 



0.0024 











1.25 











173 



Ye{Ao) = 7^\ 0.1056 
* 1.777 



The conditions for the other parameters are all given at the unification 
scale. 

5.2. Comparison of expected and measured values. We use the mod- 
ified maximal mixing boundary conditions to run the parameters and com- 
pare to the physical boundary conditions at low energy. From this analysis, 
we get that the measured Yukawa parameters for quarks are 



ld,measured(Ao) = ^ | 0.128 



0.0121 
0.128 
4.032 



and 



^ / 0.0032 

l^,measured(Ao) = TJV^ \ 0.9223 

\ 248 
For the charged leptons, we get the mass matrix 

^ , 0.000699 

>;,measured(Ao) = ^ ( « 0-147 



2.51 



Comparing these to the expected values at low energies, we see that the 
order of magnitude and form of the matrices agree. While the agreement is 
not exact, it seems that this is the closest we can get while maintaining the 
geometric constraints of the model. In order to make the agreement more 
exact, further fine tuning is required. 
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6. Conclusions 

In this paper we investigate the RGE running of the coefficients a, b, c, 5, 



e of (2.6), which appear in the asymptotic expansion of the spectral action 
functional of the noncommutative geometry model of particle physics of \TTl . 
The equations used in the renormalization group analysis are based on the 
beta function calculation of [1], for the extension of the Standard Model 
that includes right handed neutrinos with Majorana mass terms. 

We showed that the running is very sensitive to the fine tuning of the 
initial conditions at unification energy. We exhibited, as significant exam- 
ples, the different running for the default boundary conditions of [1] and the 
maximal mixing conditions, and we also showed the effect on the running 
of the coefficients of changing a single parameter in the initial conditions at 
unification. 

We then showed that a choice of boundary conditions based on the maxi- 
mal mixing, with a modified Majorana mass matrix and Higgs parameter at 
unification, satisfies all the geometric constraints on the model described in 
|13j . while at the same time gives rise to low energy values that are, within 
order of magnitude, in agreement with the expected physical values. 

We consider here the asymptotic expansion of the spectral action in the 
range of energies from the unification scale down to the electroweak scale. 
Within this range of energies, replacing the non-perturbative form of the 
spectral action with its asymptotic expansion is justified, since the error 
term is at worse of the order of A~^. However, it is known that interesting 
nonperturbative effects do arise in the spectral action, as shown in the re- 
cent results of |12j, for example in the form of a slow-roll infiation potential. 
Cosmological implications of these effects are discussed in [2^ . In terms of 
the RGE analysis considered here, we find that with our choice of modified 
maximal mixing conditions at unification, one obtains low energy values that 
are in agreement with the physical data within order of magnitude, which is 
not yet as good an agreement as one could hope for. This may be an indica- 
tion that further fine tuning of the initial conditions may achieve a better fit 
of the low energy data, or else that nonperturbative effects may play a role. 
This is not completely unlikely, considering that the nonperturbative effects 
identified in [T2j essentially appear in the coupling of Higgs and gravity and 
this in turn can affect the form of the RGE running, as observed in [T], [8], 
|32j . These questions will require further investigation. 
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